Anderson localization as position-dependent diffusion in disordered waveguides 
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We show that the recently developed self- consistent theory of Anderson localization with a 
position-dependent diffusion coefficient is in quantitative agreement with the supersymmetry ap- 
proach up to terms of the order of 1/ Qq (with qq the dimensionless conductance in the absence 
of interference effects) and with large-scale ab-initio simulations of the classical wave transport in 
disordered waveguides, at least for ^o ^ 0.5. In the latter case, agreement is found even in the 
presence of absorption. Our numerical results confirm that in open disordered media, the onset of 
Anderson localization can be viewed as position-dependent diffusion. 

PACS numbers: 42.25.Dd, 72.15.Rn 



I. INTRODUCTION 

Anderson localization is a paradigm in condensed mat- 
ter physics^. It consists in a blockade of the diffusive 
electronic transport in disordered metals due to interfer- 
ences of multiply scattered de Broglie waves at low tem- 
peratures and at a sufficiently strong disorder. This phe- 
nomenon is not unique to electrons but can manifest itself 
for any wave in the presence of disorder, in particular for 
classical waves, such as light and sound^, and, as shown 
more recently, for matter waves^. Although the absence 
of decoherence and interactions^ for classical waves is ap- 
pealing in the context of the original idea of Anderson, se- 
rious complications appear due to absorption of a part of 
the wave energy by the disordered medium^. Extracting 
clear signatures of Anderson localization from experimen- 
tal signals that are strongly affected by — often a poorly 
controlled — absorption was the key to success in recent 
experiments with microwaves^, light^ and ultrasound^. 

Classical waves offer a unique possibility of perform- 
ing angle-, space-, time- or frequency-resolved measure- 
ments with excellent resolution, the possibility that was 
not available in the realm of electronic transport. In 
a wider perspective, they also allow a controlled study 
of the interplay between disorder and interactions, as 
illustrated by the recent work on disordered photonic 
lattices^^. Interpretation of measurements requires a the- 
ory that would be able to describe not only the genuine 
interferences taking place in the bulk of a large sample 
but also the modification of these interferences in a sam- 
ple of particular shape, of finite size, and with some pre- 
cise conditions at the boundaries. Such a theory has been 
recently developed^ ^^^ based on the self-consistent (SC) 
theory of Vollhardt and Wolfle^^. The new ingredient 
is the position dependence of the renormalized diffusion 
coefficient D(y) that accounts for a stronger impact of 
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interference effects in the bulk of the disordered sam- 
ple as compared to the regions adjacent to boundaries. 
This position dependence is crucial in open disordered 
media^^. D(v) also appears in the supersymmetry ap- 
proach to wave transport ^^, which confirms that this con- 
cept goes beyond a particular technique (diagrammatic 
or supersymmetry methods) used in the calculations. 

The SC theory with a position-dependent diffusion co- 
efficient was successfully applied to analyze microwave^^ 
and ultrasonic^ experiments. The predictions of the 
theoryi^ are also in qualitative agreement with optical ex- 
periments of Storzer et al.^ . However, it remains unclear 
whether the position dependence of D is just a (useful) 
mathematical concept or if it is a genuine physical reality. 
In addition, the extent to which predictions of SC theory 
are quantitatively correct is not known. Obviously, the 
last issue is particularly important once comparison with 
experiments is attempted. 

In the present paper we compare the predictions of SC 
theory of localization with the known results obtained 
previously using the supersymmetry method^^ and with 
the results of extensive ab-initio numerical simulations 
of wave transport in two-dimensional (2D) disordered 
waveguides. We demonstrate, first, that the position- 
dependent diffusion is a physical reality and, second, that 
SC theory agrees with the supersymmetry approach up to 
terms of the order of 1/ Qq (with with qq the dimension- 
less conductance in the absence of interference effects) 
and with numerical simulation at least for ^o ^ 0-5. In 
the latter case, the agreement is found even in the pres- 
ence of absorption. 



II. SELF-CONSISTENT THEORY OF 
LOCALIZATION 

We consider a scalar, monochromatic wave u{Y)e~^^'^ 
propagating in a 2D volume-disordered waveguide of 
width w and length L ^ w. The wave field u(y) obeys 
the 2D Helmholtz equation: 



{V^^e[l^iea^3e(v)]]u(v) = {). 



(1) 



Here /c = cj/c is the wavenumber, c is the speed of the 
wave in the free space, e^ is the imaginary part of the 
dielectric constant accounting for the (spatiahy uniform) 
absorption in the medium, and (5e(r) is the randomly fluc- 
tuating part of the dielectric constant. Assuming that 
Je(r) is a Gaussian random field with a short correla- 
tion length, it is easy to show that the disorder- averaged 
Green's function of Eq. ([1]), (G(r, r')), decays exponen- 
tially with the distance |r— r'|^. The characteristic length 
of this decay defines the mean free path L In this paper 
we consider quasi- ID waveguides defined by the condi- 
tion w ^ £ <^ L. The intensity Green's function of Eq. 
(dD, C(r,rO = (47r/c)(|G(r,rO|^), obeys self-consistent 
equations that can be derived following the approach of 
Ref. 14. In a quasi-lD waveguide, all position-dependent 
quantities become functions of the longitudinal coordi- 
nate z only and the stationary SC equations can be writ- 
ten in a dimensionless form: 
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c{CC) = s{C-C). (2) 

(3) 



Here C(C, C) = {wDo/L)C{r, r'), Dq = ci/2 is the Boltz- 
mann diffusion coefficient, ( = z/L is the dimensionless 
coordinate, d{() — D{z)/ Dq is the normalized position- 
dependent diffusion coefficient, p = L/La is the absorp- 
tion coefficient (with La = ^Jit(J2 and ia = 1/^^^ 
the macro- and microscopic absorption lengths, respec- 
tively), and ^0 = ij^ /2)Ni./ L with N = kw/ir the num- 
ber of the transverse modes in the waveguide. These 
equations should be solved with the following boundary 
conditions: 



C{C,C)Tjd{C)-^C{C,C) 







(4) 



at C = and C = 1- Similarly to the 3D case^^, these con- 
ditions follow from the requirement of vanishing incoming 
diffuse ffux at the open boundaries of the sample, zq is 
the so-called extrapolation length equal to (7r/4)-^ in the 
absence of internal reffections at the sample surfacesi^. 
We will use zq = (7r/4)^ throughout this paper. When 
Eqs. (jSHI]) are solved in the diffuse regime ^o ^ 1, the 
dimensionless conductance of the waveguide is found to 
be go = {7t/2)N£/{L + 2zo)^^'^^ which is close to go for 
zo <CL. 

In the absence of absorption (/3 = 0) we can simplify 
Eq. (^ by introducing r = F(C) = J^ dCi/d{Ci): 



with the boundary conditions dU becoming 



(5) 



(6) 



and r' = F{(^)^ tq = zo/L. Equations (J5]) and (|6]) are 
readily solved: 

^/ ^/x ^ {r< + ro)(rmax + tq - r>) 

Tmax + 2ro 

where r< = min(r, t')^t^ = max(r, r') and Tmax = ^(1)- 
We now substitute this solution into Eq. (J3j) to obtain 



1 _ dr _ 



d{r) dC go 



2 ^ (r + rQ)(rmax + rQ-r) 



2ro 



This differential equation can be integrated to find r as 
a function of (. Using d{Q = {dr/dQ~^ we finally find 

d{C) = {9oVpcosh{y/pC/go) 

- [go + ro(l - p)] sinh(v^C/^o)}^ 
X {p[{~go^rof-r^p]r\ (9) 

where p is the solution of a transcendental equation 
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Solving the last equation numerically and substituting 
the result into Eq. (J9j) we can find the profile d{() at any 
go and tq = zo/L. In contrast, for /3 > Eqs. (|2]-|4]) do not 
admit analytic solution and we solve them by iteration: 
we start with D{z) = Do^ solve Eq. (|2j) numerically with 
the boundary conditions (|4]) and then find the new D{z) 
from Eq. (J3]). This procedure is then repeated until it 
converges to a solution. In typical cases considered in this 
paper the convergence is achieved after 10-20 iterations. 
The simplest object that Eqs. ([7]-[9]) allows us to study 
is the average conductance of the waveguide (g). Indeed, 
the average transmission coefficient of the waveguide is 
found as 



T = -D{L) 
1 
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where ti = F{i./ L). For the waveguide we have {g) ex T. 
A ratio that emphasizes the impact of localization effects 
is {g) /go = T/To^ where Tq is the average transmission 
coefficient found in the absence of localization effects (i.e., 
for d = 1): To = (£ + zo)/w{L + 2zo). We find 
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Simple analytic results follow for zq = 0, when go = go- 
Equation ([9j) yields 
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(13) 



and we find 



III. NUMERICAL MODEL 
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In the weak localization regime ^o ^ 1 the solution p 
of Eq. ([To]) can be found as a series expansion in powers 
of 1/go: p = 2go^ 1/3 + 2/45^o - 17/540^^ ^ ^ ^ ^^ jf 
we keep only the first term p = 2go, substitute it into 
Eq. (p!3]) and expand in powers of 1/go <C 1, we obtain 
D{z) c^Do[l- {2/go){z/L){l - z/L)]. Keeping terms up 
to 1/go in the expression for p and substituting it into 
Eqs. (dl]) and ([12]), expanding the result in powers of 1/go 
and then taking the limit of L/i ^ oo, we obtain 
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This result coincides exactly with Eq. (6.26) of Ref. [iSJ 
obtained by Mirlin using supersymmetry approach, ex- 
cept for a factor of 2 due to two independent spin states 
of electrons in Ref. [18|. We therefore proved the exact 
equivalence between SC theory and the supersymmetry 
approach for the calculation of the average conductance 
(g) up to terms of the order of l/go- 

Deep in the localized regime go <^ 1 and Eq. ([TO]) can 
be solved approximately to yield p = 1+4 exp(— l/^o) (al- 
ways for zo = and hence for go = go)- If we substitute 
this p into Eq. ([T3|) . we obtain D{z) c^ Do{ex.-p{—z/£^) + 
exp[— (L — 2;)/^]}^, where ^ = goL is the localization 
length. Equations ([M]) and ([T2]) then yield 



(16) 



where we made use of the fact that L/i ^ 1 and A^ ::^ 1. 
In contrast to Eq. ([T5|) , this result differs from the one ob- 
tained using the supersymmetry approach [see Eq. (6.29) 
of Ref. 18] . Even though the exponential decay of con- 
ductance with 1/go = L/^ — expected in the localized 
regime — is reproduced correctly, both the rate of this 
decay and the pre-exponential factor are different. We 
thus conclude that SC theory does not provide quantita- 
tively correct description of stationary wave transport in 
disordered waveguides in the localized regime. 

It is worthwhile to note that the breakdown of SC 
theory for go <^ 1 is not surprising and could be ex- 
pected from previous results. Indeed, it has already 
been noted that for the time-dependent transmission, SC 
theory does not apply after the Heisenberg time tn^'^- 
The stationary transmission coefficient T of Eq. ([TT]) 
is an integral of the time-dependent transmission T(t): 
T = Jq dt T(t)^ with the peak of T{t) around the Thou- 
less time to = I/^/tt^I^g^^- When go ~ tn/tD > 1, the 
integral is dominated hy t < tn where SC theory applies. 
The integration thus yields the correct T. However, when 
go <^ 1, tn is smaller than to and the main part of pulse 
energy arrives ai t > tn- Such long times are beyond the 
reach of SC theory, hence its breakdown for small go- 



To test the predictions of the SC model discussed in 
the previous section we solve Eq. ([1]) numerically using 
the method of transfer matrices defined in the basis of 
the transverse modes of the empty waveguide^ ^'^^. To 
this end, we represent Je(r) as a collection of M ran- 
domly positioned "screens" perpendicular to the axis z 
of the waveguide and characterized by random functions 



^N 



fv{y) = En=l Xn{y)xn{yu)- 
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5e{r) = ay^^5{z 



Zu)fu{y)- 



(17) 



Here Xn{y) = {2 / w)^ ^ '^ sin{7r ny/w) are the transverse 
modes of the waveguide and y^ are chosen at random 
within the interval {O^w). Zy represent random posi- 
tions of the screens, whereas a measures their scattering 
strength. Absorption can be included in the model by 
making ol complex. 

In the limit N -^ 00, /^^(y) becomes a delta- function 
5{y — yj,)^ mimicking a point-like scatterer. By the 
choice of fv{y) in Eq. ([TT)) we narrowed the basis to 
N right- and N left-propagating modes with real val- 
ues of the longitudinal component of the wave vector. 
Such modes are often termed "open channels" in the 
literature^^. Hence, the total transfer matrix of the sys- 
tem is a product of M pairs of 2N x 27V scattering matri- 
ces corresponding to the random screens positioned at Zj^ 
and the free space in between them, respectively^^. Be- 
cause the numerical computation of products of a large 
number of transfer matrices (~ 10^-10^ for the results 
in this paper) is intrinsically unstable, we implement a 
self-embedding procedure^^ which limits the errors in flux 
conservation to less than 10~^^ in all cases. The system is 
excited by illuminating the waveguide with N unit fluxes 
(one in each right propagating mode) and the wave field 
u(v) is computed^^^^^ for a given realization of disorder 
[see the inset of Fig. Hj^a)]. To compute statistical av- 
erages, ensembles of no fewer than 10^ realizations are 
used. 

To estimate the mean free path i of waves in our model 
system we perform a set of simulations for different disor- 
der strengths and waveguide lengths, exploring both the 
regime of classical diffusion {go > 1) and that of Ander- 
son localization (^0 < !)• The results of the simulations 
are used to compute the dimensionless conductance ^, 
equal to the sum of all outgoing fluxes at the right end 
of the waveguide, and then to study its average value 
(g) and variance var(^)^^. The dependencies of (g) and 
var(^) on go are fitted by the analytic expressions ob- 
tained by Mirlin^^ using the supersymmetry approach, 
with £ as the only fit parameter (Fig. [1])^^. The best fit 
is obtained with £ = (15.7 ± 0.2)A. In Fig. Hj^a) we also 
show Eq. ([T2]) following from SC theory. As could be 
expected from the discussion in the previous section, the 
prediction of SC theory coincides with both the results of 




FIG. 1: The average (a) and the variance (b) of the conductance g of disordered waveguides supporting A^ = 10 (circles) and 
A/" = 20 (squares) modes are shown versus the inverse of go. The sohd hues marked as SUSY are fits using Eq. (6.23) of Ref. 
\m , derived using the supersymmetry approach, with -^ = 15. 7 A as the only fit parameter. The soUd line marked as SC in (a) is 
obtained using the self- consistent theory [Eq. ()12p ]. Inset in (a): for a given realization of disorder, wave "trajectories" found 
by connecting local Poynting vectors are superimposed on the distribution of intensity |i^(r)|^ in a disordered waveguide with 
w = 10.25A and L = 50A. Only trajectories that traverse the waveguide are shown. 



the supersymmetry approach and numerical simulations 
only for large go > 0.5. 



IV. POSITION-DEPENDENT DIFFUSION 
COEFFICIENT 

The wave field u{r) that we obtain as an outcome of 
the numerical algorithm allows us to calculate the energy 
density >V(r) and flux J(r)^^: 

—kc Im [w(r) Vu(r)] . 



vy(r) 

J(r) 



\Vu{r)f 



(18) 
(19) 



These two quantities formally define the diffusion coeffi- 
cient D{z) which, in general, may be position-dependent: 



D{z) 



(Mr)) 



(20) 



where the averages (...) are taken over a statistical en- 
semble of disorder realizations as well as over the crossec- 
tion of the waveguide. Eq. ([2Q|) can be used only at dis- 
tances beyond one mean free path i from the boundaries 
of the random medium because more subtle propagation 
effects of non-diffusive nature start to be important in 
the immediate vicinity of the boundaries^. 

We first consider non-absorbing disordered waveguides 
described by e^ = in Eq. ([1]) and real a in Eq. (p!7|) . 
In Fig. [2] we compare numerical results for D{z) with 
the outcome of SC theory for waveguides of different 
lengths but with statistically equivalent disorder. Quan- 
titative agreement is observed for L = 100-800A, corre- 
sponding go ^ 0.3-2. For the longest of our waveguides 



{L = 1600 A, go ^ 0.16), deviations of numerical results 
from SC theory start to become visible in the middle of 
the waveguide, which is particularly apparent in the log- 
arithmic plot of Fig. [2fb). The mean free path i = 17. 5 A 
corresponding to the best fit of SC theory to numerical 
results is only about 10% higher than £ = 15. 7A obtained 
from the fits in Fig. [H 

We checked that the results of numerical simulations 
are not sensitive to the microscopic details of disor- 
der: D{z) obtained in two runs with different scattering 
strengths a and different scatterer densities, but equal 
mean free paths i turned out to be the same. 



V. EFFECT OF ABSORPTION 

The linear absorption is modeled by introducing a non- 
zero Ca in Eq. ([!]) and making a in Eq. (fT7|) complex. 
A link between e^ and a can be established using the 
condition of flux continuity. Indeed, for continuous waves 
considered in this work the continuity of the flux leads to 



(V-J(r)) = (c/£a)(W(r)), 



(21) 



where ia — l/kca. We checked that within numerical ac- 
curacy of our simulations the proportionality factor c/£a 
indeed remains constant independent of z. Therefore, 
Eq. (|2T]) allows us to determine the microscopic absorp- 
tion length ia as c(W(r))/(V- J(r)) obtained numerically 
at a given a. 

Figure [3] demonstrates the effect of absorption on the 
position-dependent diffusion coefficient for a waveguide 
of length L = 400A, which is about 25 mean free paths. 
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FIG. 2: (a) Position-dependent diffusion coefficient D[^z) in 
2D waveguides supporting the same number N — 10 of trans- 
verse modes (width w — 5. 25 A) but having different lengths 
L. Disorder is the same for ah lengths. Symbols show the 
results of numerical simulations, whereas solid lines are ob- 
tained from the self- consistent theory with the mean free 
path i — 17.5 A. Dashed lines show the approximate results 
for ^0^1 (shown for L — 100 A) and ^o ^ 1 (shown for 
L = 1600A), with D(0) substituted for D^, see text, (b) 
Same as (a) but in the logarithmic scale. 



For this waveguide ^o — 1-3 and the localization cor- 
rections are important. We observe that absorption 
suppresses the localization correction to the position- 
dependent diffusion coefficient. This clearly demon- 
strates that the absorption nontrivially affects the trans- 
port by changing the way the waves interfere. Never- 
theless, we observe good agreement between numerical 
results (symbols) and SC theory (solid lines). The pre- 
dictions of SC theory start to deviate from numerical 
results only for strong absorption {La/L < 0.4). Once 
again, the mean free path I = 17.1 A obtained from the 
fit of SC theory to the lower curve of Fig.[3]is within 10% 
of the value estimated from the variance of dimensionless 
conductance. 
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FIG. 3: The effect of absorption on the position-dependent 
diffusion coefficient. Symbols are results of numerical simula- 
tions in a 2D waveguide of length L = 400A, width w = 10.25A 
{N = 20) and several values of the macroscopic absorption 
length La indicated on the figure. Lines are obtained from 
SC theory with i = 17.1 A adjusted to obtain the best fit for 
the case of no absorption (lower curve). Dashed line shows 
D(z) following from the self- consistent theory with the same 
i = 17. 5A as in Fig. [2] and illustrates the sensitivity of D{z) 
to the exact value of i. 



VI. CONCLUSIONS 

Two important results were obtained in this work. 
First, we convincingly demonstrated that the position- 
dependent diffusion coefficient is not an abstract mathe- 
matical concept but is a physical reality. The results of 
numerical simulations of scalar wave transport in disor- 
dered 2D waveguides unambiguously show that the on- 
set of Anderson localization manifests itself as position- 
dependent diffusion. The reduction of the diffusion co- 
efficient D{r) is much more important in the middle of 
an open sample than close to its boundaries, in agree- 
ment with predictions of the self-consistent theory of lo- 
calization. Second, we established that for monochro- 
matic waves in 2D disordered waveguides predictions of 
the self-consistent theory of localization are quantitatively 
correct provided that the dimensionless conductance in 
the absence of interference effects go is at least larger 
than 0.5. Moreover, the self-consistent theory yields a 
series expansion of the average conductance (g) in pow- 
ers of 1/go that coincides exactly with the expansion ob- 
tained using the supersymmetry method^^ up to terms 
of the order of l/g^. This was not obvious a priori be- 
cause of the numerous approximations involved in the 
derivation of self-consistent equations^^. The agreement 
between theory and numerical simulations is good in the 
presence of absorption as well, which has a particular im- 
portance in the context of the recent quest for Anderson 
localization of classical waves that heavily relies on con- 
frontation of experimental results with the self-consistent 
theory "ii^ii^. Deep in the localized regime (^o < 0-5), 
the self-consistent theory loses its quantitative accuracy, 
but still yields qualitatively correct results (exponential 



decay of conductance with the length of the waveguide 
and of the diffusion coefficient D with the distance from 
waveguide boundaries). It would be extremely inter- 
esting to see if the ability of the self-consistent theory 
to provide quantitative predictions still holds in three- 
dimensional systems where a mobility edge exists. In 
particular, the immediate proximity of the mobility edge 
is of special interest. 

Note added. After this paper was submitted for publi- 
cation, a related preprint appeared^^. In particular, the 
authors of that work show that the self-consistent theory 
does not apply to ID disordered media, which is consis- 
tent with our results because qq ~ (./ L is always small in 
ID, provided that the condition L ^ I assumed in this 



paper is fulfilled. 
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